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Statement  of  the  problem  studied 


This  research  project  focussed  on  the  enhancement  of  current  vehicle  sim¬ 
ulation  capabilities  of  TARDEC.  Problems  addressed  included  development 
of  new  numerical  algorithms  for  efficiently  simulating  mechanical  systems 
with  low-amplitude,  high-frequency  vibrations,  and  development  of  numeri¬ 
cal  methods  and  models  for  impacts  and  collisions. 


2  Summary  of  the  most  important  results 


High-frequency  oscillations  in  nonlinear  ODE/DAE  systems  are  a  problem 
in  vehicle  simulation  because  following  the  oscillations  necessitates  the  use 
of  an  extremely  small  timestep.  However,  many  of  the  oscillations  are  not 
important  for  the  overall  numerical  solution.  We  have  shown  that  methods 
based  on  local  linearization  can  fail  because  the  local  eigenstructure  of  the 
problems  oscillates  at  the  high  frequency.  Experiments  have  demonstrated 
that  certain  implicit  methods  combined  with  automatic  stepsize  control  can 
damp  out  the  oscillation  safely,  in  regions  where  its  amplitude  is  too  small 
to  be  important.  The  usual  stepsize  selection  strategies  must  be  modified 
so  that  they  are  correct  for  the  limiting  high-index  DAE.  There  is  also  some 
theory  to  support  this  technique  for  the  equations  of  motion.  However,  once 
the  stepsize  is  increased,  problems  with  Newton  iteration  convergence  again 
restrict  the  timestep.  These  problems  are  also  due  to  the  rapidly  chang¬ 
ing  local  eigenstructure.  Some  formulations  of  the  equations  of  motion  are 
more  advantageous  than  others  in  terms  of  Newton  iteration  Convergence  for 
highly  oscillatory  systems.  A  coordinate-split  (CS)  method  has  been  devel¬ 
oped  that,  together  with  a  modified  Newton  (CM)  iteration  is  particularly 
effective.  Numerical  results  for  a  number  of  highly  oscillatory  multibody  sys¬ 
tems  demonstrate  that  the  new  method  is  particularly  effective  for  highly 
oscillatory  systems  where  the  oscillation  is  small  and  can  be  damped.  Re¬ 
cently  developed  theory  explains  the  Newton  convergence  results.  Investiga¬ 
tion  of  the  high-frequency  oscillation  problem  was  motivated  by  discussions 
with  Roger  Wehage  and  Jim  Over  holt  (TARDEC). 
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Modeling  impact  of  bodies  or  particles  during  a  physical  system  simulation 
is  problematic  and  may  cause  many  difficulties  in  the  numerical  solution.  On 
the  other  hand,  it  can  better  define  the  dynamics  of  some  systems,  such  as  in 
track  vehicle  simulation.  It  is  well-known  that  a  purely  algebraic  treatment 
of  the  rigid  impact  problem  can  lead  to  incorrect  results  in  some  cases.  To 
remedy  this,  Keller  proposed  a  set  of  evolution  equations  which  can  be  used 
for  non-trivial  three-dimensional  impact  problems.  For  the  hypothesis  of 
impact,  Stronge  proposed  a  resititution  model  that  overcomes  the  difficulty 
of  using  either  Newton’s  or  Poisson’s  stopping  criterion.  A  rigid  impact 
model  between  the  road-wheel  and  track  has  been  developed  using  Keller’s 
evolution  equations.  The  coefficients  of  friction  and  restitution  are  the  only 
input  parameters.  By  adjusting  the  coefficients,  accurate  tensional  forces 
applied  to  the  track  can  be  computed  using  the  tangential  impulse.  Two- 
body  impact  evolution  equations  were  generalized  for  the  constrained  multi- 
rigid- body  impact  model,  e.g.,  rigid  impact  between  the  road-arm  and  road- 
wheel  composite  body  and  the  chain  of  track  segments.  Compared  to  the 
conventional  track  model,  where  one  has  to  guess  the  soil  deformation,  and 
then  compute  the  shear  and  normal  forces  between  the  ground  and  track,  this 
approach  yields  a  more  accurate  and  effective  computational  scheme  to  treat 
the  complex  track  model.  This  work  has  been  in  collaboration  with  G.P. 
Mac  Sithigh  at  the  Mechanical  and  Aerospace  Engineering  and  Engineering 
Mechanics  department.  University  of  Missouri-Rolla. 
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